Hexagonal pnictide SrPtAs: superconductivity with locally broken inversion 
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Unlike other pnictides, SrPtAs has a hexagonal structure, containing layers with As-Pt atoms 
that form a honeycomb lattice. These layers lack inversion symmetry which allows for a spin-orbit 
coupling that we show has a dramatic effect on superconductivity in this material. In particular, 
for conventional s-wave superconductivity in SrPtAs, both the spin susceptibility and the param- 
agnetic limiting field are enhanced significantly with respect to that usually expected for s-wave 
superconductors. SrPtAs provides a prime example of a superconductor with locally broken inver- 
sion symmetry. 



The superconducting pnictides [l[ present a fasci- 
nating class of materials that highlight the interplay 
between electronic correlations, superconductivity, and 
magnetism in a multi-orbital system |2[ . SrPtAs is a new 
member to this family with a unique feature: the As-Pt 
atoms in a single layer form a honeycomb lattice Q, see 
Fig. [TJ It is natural to ask if the new lattice structure 
can have a consequence on superconductivity. Here we 
argue that it does. In particular, even though SrPtAs 
has a center of inversion symmetry, the broken inver- 
sion symmetry inherent to a single As-Pt layer has non- 
trivial consequences. Assuming that SrPtAs is a spin- 
singlet superconductor, we show that it is expected to 
have a non- vanishing spin susceptibility at zero tempera- 
ture with a magnitude that is a significant portion of the 
normal state spin-susceptibility. We further show that it 
is likely to have a critical field larger than the paramag- 
netic limiting field. 

Of particular importance to this work is that a single 
As-Pt layer lacks a center of inversion symmetry. The 
unit cell of SrPtAs contains two inequivalent As-Pt lay- 
ers that are related by inversion symmetry, see Fig. [T] 
As shown below, SrPtAs also has a spin-orbit coupling 
(SOC) that is larger than the inter-layer coupling. This 
combination of broken local inversion and strong spin- 
orbit coupling allows SrPtAs to be a good candidate for 
superconductivity with local inversion-symmetry break- 
ing. We use this term to refer to the fact that physical 
properties usually associated with non-centrosymmetric 
superconductivity appear in SrPtAs, despite the pres- 
ence of a center of inversion symmetry. For spin-singlet 
superconductors, these properties include an enhanced 
paramagnetic depairing field and a non-vanishing spin- 
susceptibility at zero temperature 0-0] • In the following, 
we initially present the electronic structure of SrPtAs and 
then turn to an examination of the superconducting state 
in this material. 



First-principles calculations were performed using the 
highly precise full-potential linearized augmented plane 
wave (FLAPW) method We have used the exper- 

imental lattice constants a = 4.24A and c = 8.98A(Tl| 
and a cutoff of 186 eV for basis functions. The local 
density approximation (LDA) is used for the exchange- 
correlation as parametrized by Hedin and Lundqvist[12j, 
and spin-orbit coupling has been calculated using a 
second- variational treatment [l3| . Fig. [2] and Fig. [3] show 
the results of LDA calculations with and without spin- 
orbit coupling. Energy bands near the Fermi level orig- 
inate from Pt 5d and As Ap orbitals. Specifically, the 
Fermi surface sheets labeled a and b in Fig. [3] (a) stem 
from Pt d xy , d x 2_ y 2, As p x , and p y orbitals while that 
labeled c stems from Pt d XZl d yz and As p z orbitals. Our 
results without spin-orbit coupling agree with those of 
Ref. 14|. Note the qualitative changes when spin-orbit 
coupling is added. In particular, the spin-orbit coupling 
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FIG. 1. (a) Structure and (b) Brillouin zone in SrPtAs. Red, 
blue, and grey spheres denote Pt, As, and Sr atoms, respec- 
tively. 
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FIG. 2. Energy bands of SrPtAs (a) without and (b) with 
spin-orbit coupling. Zero energy represents the Fermi level. 
Indices a, b, and c in (a) represent the three bands crossing 
the Fermi level. 

leads to appreciable changes in the band structure along 
the symmetry lines of H — A and L — H . Also of rele- 
vance is the difference between the bands along the sym- 
metry lines H — A — L and K — Y — M when there is 
no spin-orbit coupling (Fig. [2]). This difference is due to 
inter-layer coupling between the As-Pt layers. This cou- 
pling vanishes for symmetry reasons in the plane given 
by k z = 7r/c. The band structure reveals that the band 
splittings due to spin-orbit coupling are comparable to or 
larger than those due to inter-layer coupling. A fact that 
plays an important role in the superconducting state. 

To understand the bands stemming from the LDA cal- 
culations, it is useful to consider initially a single As-Pt 
layer. A key point is that this layer does not have a cen- 
ter of inversion and, therefore, a spin-orbit coupling term 
of the form 

fe,s,s' 

exists, where c, (c^ .) creates (annihilates) an elec- 
tron with momentum k and pseudo-spin s in layer i, 
a denote the Pauli matrices, and on is the layer i 
spin-orbit coupling energy. Time-reversal symmetry im- 
poses g k — —g_ k - The invariance of the Hamilto- 
nian under point-group operations leads to the require- 
ment g k = G s g G _ 1 ^,, where G s is the rotation matrix 
for a pseudo-vector and G is the rotation matrix for 
a vector. The point group of a As-Pt layer is D^. 
Invariance under the mirror symmetries with normals 
along the c axis and along the Pt-Pt bond imply that 
g(k x ,0,0) = g(fc x ,0,7r/e) = (here a Pt-Pt bond is 
taken to be along the y axis) . This reveals itself for bands 
along the A to L direction in the Brillouin zone, where 
there is no spin-orbit splitting, see Fig. 3(b). Within a 
tight-binding approach, we find g k — i^ i sin(fe • Tj), 
where Tj are the translation vectors T\ = (0, a, 0), 
T 2 = (V3o/2, -a/2,0), and T 3 = (-■ \/3a/2, -a/2, 0). 
This form of spin-orbit coupling can be found for all 
bands stemming from Pt d orbitals by including hop- 
ping to neighboring As p orbitals and by including on-site 
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FIG. 3. Cross sections of the Fermi surface of SrPtAs with and 
without spin-orbit coupling. Figure (a) [(c)] is for k z — tt/c 
and no spin-orbit [with spin-orbit], while Figure (b) [(d)] is for 
k z — and no spin-orbit [with spin-orbit]. Indices a, b, and c 
in (a) represent the three bands crossing the Fermi level. 

L • S for both As and Pt sites. We note that this spin- 
orbit coupling is similar to that used by Kane and Mele 
to discuss the quantum spin hall effect in graphene fljj . 
Symmetry also allows for g x and g y to be non-zero. How- 
ever these components must be odd in k z and, within a 
tight-binding analysis, are only found by including hop- 
ping along the c-axis. Given the much weaker dispersion 
of the bands along k z relative to the in-plane dispersion, 
we expect that g x and g y are much smaller than g z and 
we will only include g z in the following. 

The analysis in the previous paragraph applies to a 
single As-Pt layer. The two inequivalent As-Pt layers 
are related by inversion symmetry, consequently, oti is 
of opposite sign for the two layers, i.e., a; = (— l) l a. To 
complete the description for the solid, a coupling between 
the two inequivalent layers is required. We take this to 
be e c (fc) (symmetry requires this to vanish for k z — tt/c). 
Provided there are no band degeneracies other than spin 
and layer degeneracies, a generic Hamiltonian for SrPtAs 
is then 

Ho =^# t (fe)|[ei(fc) - /i]cr T + a(k)a z r z 
k 

+ Re[e c (k)]a T x + Jm[e c (fe)]cr r H |*(fc), (2) 

where *(fe) = (c kn , c kil , c fet2 , c ki2 ) T , Oi (n) are Pauli 
matrices that operate on the pseudo-spin (layer) space 
and a(k) — ag z (k). This Hamiltonian can be diagonal- 
ized with resulting dispersion relations e±(k) = ei(k) ± 
^/|e c (k)| 2 + a 2 (k) and each state has a 2-fold Kramers 
degeneracy. To gain an intuition for the terms ap- 
pearing in this Hamiltonian, we state the results for a 
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simple tight-binding theory (note that below we keep 
these terms arbitrary). This yields ei(k) = ii(cosfc • 
Tx + cosfc • T 2 + cosfc • T 3 ) + t c2 cos(ck z ) and e c (k) = 
t c cos(fc 2 c/2)(l + e" ik T3 +e ik T2 ) (note that g(k) and Ti 
are given above). 

From the point of view of superconductivity, the strong 
spin-orbit coupling is of interest. In the limit that the 
inter-layer coupling vanishes, we have two uncoupled 
non-centrosymmetric systems. It is known that in non- 
centrosymmetric spin-singlet superconductors the spin- 
susceptibility and the paramagnetic limiting field are sig- 
nificantly enhanced, if the spin-orbit coupling strength is 
much larger than the superconducting gap [^-[Sj] - Given 
the large spin-orbit coupling relative to the inter-layer 
coupling, it is conceivable that the behavior of supercon- 
ducting SrPtAs resembles that of a non-centrosymmetric 
material. For this reason we calculate both the spin- 
susceptibility and the limiting field assuming that super- 
conductivity in SrPtAs is spin-singlet (this is a reasonable 
assumption comparing with other pnictides superconduc- 
tors [2j). To be concrete we assume intra-layer s-wave 
pairing with an interaction 



n sc = -v 



k,k' ,i,s,s' 



j 

fc.s 



-k's'i c k' 



(3) 



Note that our results do not depend qualitatively on this 
choice. 

For a system described by the Hamiltonian ([2]) and 
(J3|) , the susceptibility in the superconducting and normal 
state can be calculated using [15| 

Xij = -Ms^y y tr[cr i G(fc,cJ, l )cr 3 G(fc,cJ, l ) 
n k 

-aiF{k,u n )(7jF\k,u n )] (4) 

with G(k,u! n ) and F(k,ui n ) the normal and anomalous 
Green's functions in the Matsubara formulation. Note 
that even for this 'one-band' formulation, the Green's 
functions are 4x4 matrices, so that the trace runs 
over both, layer and spin index. In the notation of the 
Hamiltonian ([2]) there are essentially three bands cross- 
ing the Fermi energy in SrPtAs [labeled a, b, and c in 
Figs. [2ja) and EJa)] and we can generalize the above ex- 
pression to 



Xij 



(5) 



where the sum runs over the three bands v = a,b, c. 
Below we calculate the susceptibility separately for each 
band using Eq. ((4]). 

In the normal state, F v (k,uj n ) = 0, and we find for 
fields parallel to z, 



X 



dn F (eUk)) 



k,i=± k 



= Y,X° P (k,v), (6) 



where n F (e) is the Fermi distribution function as a func- 
tion of energy e and x l p(k, v) denotes a Pauli susceptibility 
for band v. For fields in plane, we find 



|e£(fc)|a + [o»(fc)]2 



with the van Vleck susceptibility 



(7) 



X°vv( k ^) 



_ 2 rM4(*0) 

z t x B\ 



n F (e»_{k)) 



}• (8) 



The Pauli susceptibility describes intra-band processes 
and at low temperatures is proportional to the density 
of states at the Fermi level. The van Vleck suscepti- 
bility describes inter-band processes. For the supercon- 
ducting states in the limit ,/|e£(fc)| 2 + [a"{k)] 2 > A" , we 
recover the expressions given in Eqs. (J6)) and ([7]). where 
for the Pauli susceptibility, we have to replace e± (fe) with 
El(k) = s/[e^(k)] 2 + [A"] 2 . The Pauli susceptibility con- 
tribution thus vanishes due to the opening of the super- 
conducting gap, while the van Vleck susceptibility is un- 
changed by superconductivity, even at T = 0. Conse- 
quently, Xz C will behave like that expected for a conven- 
tional spin-singlet superconductor while Xj3 will have a 
large spin susceptibility, even at T — 0. To demonstrate 
this, Fig. 0] shows the ratio of XxzM m the supercon- 
ducting phase at T — to the normal state in-plane 
susceptibility x\ z { v ) as a function of a v ' jt v c for the three 
bands v = a,b,c (where a u is the spin-orbit strength 
and t v c is the interlayer coupling strength). These values 
were determined using simple tight-binding calculations 
for the three Fermi surface sheets a, b, c. Estimating the 
ratios a" ' jt v c from the band structure, we find values of 
X S ° z {v)lx\ z {v) = 0.11,0.42, and 0.91 for Fermi surface 
sheet v — a, b, and c, respectively. Consequently, we 
expect that a sizable portion of the normal state suscep- 
tibility will exist in the limit T — ¥ for in-plane magnetic 
fields. We note that related behavior has recently been 
predicted for the spin-susceptibility in an examination of 
the crossover from non-centrosymmetric to centrosym- 
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metric superconductivity in multi-layer systems 

The enhanced susceptibility suggests that the Pauli 
limiting field will also be enhanced when the field in the 
basal plane. To calculate this, we include the Zeeman 
field Uz = Z}fe, s , s /,iW-BH ■ cr 8 ^,^ kai c ks ,. and orient the 
field in the basal plane. Within weak-coupling theory 
and assuming \/|e c (fc)| 2 + a(k) 2 2> gn-BH (which is well 
supported by LDA results) , we find the following expres- 
sion for T c as a function of h = <7/is|H|: 



In 



(—) 



Re 



where 



h(k) 



h\n c (k)\ 
2nk b T c 



V|e c (fc)| 2 + <* 2 (fe)- 



(10) 
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FIG. 4. Spin susceptibility at T = in the superconducting 
state normalized with respect to the normal state susceptibil- 
ity for the three bands crossing the Fermi surface as a func- 
tion of a" /t". The black squares denote approximate values 
for these three bands. 



perconductivity. In particular, the lack of an inversion 
center in the As-Pt honeycomb lattice layers, combined 
with strong spin-orbit coupling, allows a significant en- 
hancement of both the Pauli limiting field and the spin 
susceptibility for spin-singlet superconductivity. SrPtAs 
provides an ideal example of superconductivity with lo- 
cally broken inversion symmetry. 
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T c q is the transition temperature for h = 0, (f)^ means 
an average of / over the Fermi surface, Re means real 
part, and ^>(x) is the digamma function. The band index, 
v, is omitted for brevity. In the limit that e c = 0, we find 
that T c is independent of h. This is in agreement with the 
result predicted and observed for non-centrosymmetric 
superconductors in the limit of large spin-orbit coupling 
[H, Also, in the limit that a = 0, we find the usual 
expression for the Pauli limiting field. At T c = 0, using 
Eq. ©, we find that the Pauli limiting field is given by 

,1N = / h\e c (k)\ \ 

\2J \ W B T c0v /|e c (fc)| 2 + a 2 (fc) /fe 

Using tight-binding calculations, we estimate that the 
enhancement of the Pauli limiting field, hp/hpo (where 
hpo is the limiting field when a — 0), takes the values 
1.1, 1.8, and 7.4 for Fermi sheet a, b, and c, respectively. 
Provided that the orbital upper critical field is sufficiently 
large, it should be possible to observe an enhanced Pauli 
limiting field. For fields along the c-axis, a usual Pauli 
suppression is expected. 

Finally, we point out that in addition to spin-singlet 
pairing, a spin-triplet component will also appear[l^|. In 
particular, in a given layer a spin-triplet component with 
d(k) along the direction of g(k) exists, such that it has 
opposite sign in the two inequivalent layers of SrPtAs. 

In conclusion, we have shown that the unique struc- 
ture in the pnictide SrPtAs has non-trivial effects on su- 
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